基于Hellinger-Reissner变分原理的应变梯度杂交元设计

从一般的偶应力理论出发，基于Hellinger-Reissner变分原理，通过对有限元 离散体系的位移试解引入非协调位移函数，得到了偶应力理论下有限元离散系统的能量相容 条件，并由此建立了应变梯度杂交元的应力函数优化条件. 根据该优化条件，构造了一 个C0类的平面4节点梯度杂交元，数值结果表明，该单元对可压缩和不可压缩状态的 梯度材料均可给出合理的数值结果，再现材料的尺度效应.     

复合纤维增强混凝土阻尼测试装置开发与试验研究

纤维与聚合物的掺入可以明显改善混凝土材料的阻尼性能。本文首先给出了正弦交变激励下粘弹性材料三点弯曲梁阻尼特性关系,其次首次自主开发了大尺寸材料的阻尼性能测试装置,然后利用开发的装置在频率(0.5~2.0Hz)条件下测定了6种不同配比复合纤维增强阻尼混凝土的损耗因子与储存模量,最后对纤维的阻尼增强机理进行了初步探讨。试验结果表明:复合纤维增强阻尼混凝土与素混凝土相比,提高了混凝土的损耗因子80%~200%。主要原因是聚合物分子在外力作用下的内耗增加了普通混凝土的阻尼能力,而纤维的阻尼增强机理在于纤维的掺入增加了纤维与水泥基材的界面摩擦力。     

一种求解层合圆柱壳的杂交应力单元

根据修正的余能原理，推导出一种求解复合材料层合圆柱壳的杂交应力单元。取用六面体等参单元，此单元反映了各层材料性质不同及应力分布沿整个厚度不连续现象，同时计入横向剪切变形和法向挤压变形，适用于厚层壳体。文章通过实例说明此单元能准确求出各层内的应力值，实用价值高。     

结构形状优化与智能模型化

阻碍形状优化软件广泛应用的障碍之一是依据自然设计度量描述和建立设计模型、分析模型、优化模型以及实现三个模型之间的转换.本文称这一困难为结构形状优优设计软件的适用性.本文提出了一种方法用来动态地确定平面连续体结构形状优化过程中的边界,应用基于人工智能方法的启发式规则与技术,自动生成由设计单元法表示的几何模型,也就是将一个结构自动剖分成若干个大的四边形映射单元.这些大单元对于进一步的网络生成是必要的,同时也是向全自动计算机辅助形状优化系统前进的重要一步.     

计算钝感炸药爆轰过程的一种新反应率模型

根据钝感炸药爆轰过程中含有激化过程和慢反应,建立了一种新的反应率模型。与其他反应率模型相比,这种反应率模型可以应用于较粗的网格。在每厘米50个网格条件下,炸药驱动铝飞片和钽飞片的自由面计算结果与实验结果很接近。同时,应用此反应率模型计算了钝感炸药驱动LiF过程,在每厘米50个网格的条件下,炸药与LiF间速度的计算值同实验值接近,且误差随网格尺寸变小而变小。这些表明,此反应率模型能够在较粗的网格条件下,比较准确地描述钝感炸药驱动飞片过程,有利于在工程实际中应用。     

人颅骨结构的动力参数

为了深入探讨人头部受撞击时颅脑损伤的力学计算模型,本文对颅骨结构的功力参数连行了测试和有限元分析。结果表明,人颅骨结构属于小阻尼结构,进行动力分析时可以抽象成弹性阻尼薄壳结构。     

Hahn-Tsai复合材料的非线性杂交应力有限元方法

成功建立了Hahn-Tsai复合材料模型的非线性杂交应力有限元方程,采用Newton-Raphson迭代法求解结构的非线性位移方程。在迭代过程中,为了提高计算效率可采用简单迭代法由节点位移求解单元应力场。但是,当载荷增加到一定程度以后,非线性应力场由于循环迭代而无法收敛,显然,一般的加速方法不能解决这种循环迭代的发散问题。因此,本文发展了一种确实有效的非线性应力场迭代新方法,在不增加计算工作量的情况下,不仅极大地提高了收敛速度,而且对于较大载荷也能够很好地收敛,从而解决了大载荷下非线性杂交元方法失败的关键问题。数值算例表明该方法是确实可行的。     

有限元网格自动生成及自动调整技术

本文给出了一种任意平面域内三角形网格自动生成及自动调整算法，并编制了相应的程序。通过对带Ｕ形槽的三点变曲试件的分析，证明本文算法是可行的，并且具有稳定性能好、收敛速度快的优点。     

Numerical solution of Navier–Stokes equations using multiquadric radial basis function networks

A numerical method based on radial basis function networks (RBFNs) for solving steady incompressible viscous flow problems (including Boussinesq materials) is presented in this paper. The method uses a ‘universal approximator’ based on neural network methodology to represent the solutions. The method is easy to implement and does not require any kind of ‘finite element‐type’ discretization of the domain and its boundary. Instead, two sets of random points distributed throughout the domain and on the boundary are required. The first set defines the centres of the RBFNs and the second defines the collocation points. The two sets of points can be different; however, experience shows that if the two sets are the same better results are obtained. In this work the two sets are identical and hence commonly referred to as the set of centres. Planar Poiseuille, driven cavity and natural convection flows are simulated to verify the method. The numerical solutions obtained using only relatively low densities of centres are in good agreement with analytical and benchmark solutions available in the literature. With uniformly distributed centres, the method achieves Reynolds number Re = 100 000 for the Poiseuille flow (assuming that laminar flow can be maintained) using the density of , Re = 400 for the driven cavity flow with a density of and Rayleigh number Ra = 1 000 000 for the natural convection flow with a density of . Copyright © 2001 John Wiley & Sons, Ltd.     

A two-grid method for elliptic problem with boundary layers

A two-grid method for the elliptic equation with a small parameter ε multiplying the highest derivative is investigated. The difference schemes with the property of ε-uniform convergence on a uniform mesh and on Shishkin mesh are considered. In both cases, a two-grid method for resolving the difference scheme is investigated. A two-grid method has features that are concerned with a uniform convergence of a difference scheme. To increase the accuracy, the Richardson extrapolation in two-grid method is applied. Numerical results are discussed.